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ABSTRACT 

We present radio continuum observations at 7 mm made using the Very Large 
Array towards three massive star forming regions thought to be in very early 
stages of evolution selected from the sample of Sridharan et al. (2002). Emission 
was detected towards all three sources (IRAS 18470-0044, IRAS 19217+1651 and 
IRAS 23151+5912). We find that in all cases the 7 mm emission corresponds to 
thermal emission from ionized gas. The regions of ionized gas associated with 
IRAS 19217+1651 and IRAS 23151+5912 are hypercompact with diameters of 
0.009 and 0.0006 pc, and emission measures of 7.0 x 10^ and 2.3 x 10^ pc cm~^, 
respectively. 

Subject headings: nebulae: HII regions — stars: formation — stars: massive 



1. Introduction 



One of the key challenges at present in the field of star formation is to understand 
how massive stars (M > 10 Mq) form. Our present understanding of star formation is 
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primarily based on the observations of low-mass stars. In the current paradigm (Shu, Adams, 
& Lizano 1987; Shu et al. 1993) the formation of low-mass stars is characterized by an 
accretion phase, in which a central protostar and a circumstellar disk form surrounded by 
an infalling envelope of dust and gas, followed by a phase in which the protostar deposits 
linear and angular momentum, and mechanical energy into its surroundings through jets 
and molecular outflows. Although this paradigm has been very successful in explaining 
what is observationally known about the formation of low-mass stars (e.g, Lada 1991), its 
applicability to the formation of massive stars remains arguable. In particular, alternative 
mechanisms such as the merging of lower mass protostars to form a massive protostar have 
received serious consideration lately (e. g. Bonnell, Bate, & Zinnecker 1998; Portegies Zwart 
ct al. 1999), and there is evidence of strong dynamical interaction and possibly merging 
deduced from the proper motions of some young massive stars in the Orion KL/BN region 
(Gomez et al. 2005). 

Whether the massive protostars are formed by accretion or merging remains controver- 
sial. If massive O stars are formed by accretion, we expect that disks and jets will be present 
in their earliest stages of evolution, as in the case of low mass stars. Due to their large 
luminosities, hypercompact Hll regions are also expected to be formed while the protostar 
is still undergoing accretion (Keto 2002, 2003). On the other hand, if formed by random 
coalescence of lower-mass stars neither disks, jets or hypercompact HII regions are expected. 
Thus, the search for jets, disks and hypercompact HII regions toward massive YSOs is crucial 
to understand their formation process. 

The presence of jets and disks toward massive YSOs still lacks solid grounds. CoUimated 
outflows and/or disks have been found in a handful of B-type protostars: IRAS 18162-2048 
(£ ~ 1.7 X 10^ Lq; Marti, Rodriguez, & Reipurth 1993); Cepheus A HW2 (£ ~ 1 x 10^ Lq] 
Rodriguez et al. 1994); IRAS 20126+4104 (£ ~ 1.3 x 10^; Cesaroni et al. 1997); G192.16- 
3.82 (£ - 3 X 10^; Shepherd et al. 1998; Devine et al. 1999; Shepherd, Claussen, & Kurtz 
2001); AFGL 5142 {C ~ 3.0 x 10^ Lq; Zhang et al. 2002). None of these objects exceeds 
a bolometric luminosity of 2 x IQ'^Lq and are thus BO ZAMS or lower luminosity objects. 
The most luminous 0-type protostars associated with a jet and a collimated molecular 
outflow is IRAS 16547—4247, a luminous infrared source with a bolometric luminosity of 
6.2 X 10^LQ(Garay et al. 2002, 2007; Rodriguez et al. 2005). There is, however, no known 
unquestionable case of a circumstellar disk associated with a massive 0-type protostar. 
The search for disks around independent high-mass protostars is inherently difficult. The 
photoevaporation timescale of disks by the radiation from the central star (Yorke & Welz 
1996; Hollenbach 1997) can be relatively short, and disks may be destroyed before the 
dispersal of the molecular cloud core. If other OB stars in the region have been formed 
before, the object studied may be embedded in bright free-free emission that will difficult 
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the search for the relatively weak mm emission from the disk. 

We report here 7 mm continuum observations, made with the Very Large Array, towards 
three luminous (L > 6 x 10^ Lq) objects, selected from the sample of high-mass protostellar 
candidates of Sridharan at al. (2002), that are associated with powerful bipolar CO out- 
flows (indicating the presence of a coUimated outflow and thus suggesting the presence of 
a disk). The main goal of these observations was to search for either circumstellar disks or 
hypercompact regions of ionized gas around 0-type protostars. 



The observations were made with the Very Large Array (VLA) of the National Radio 



good weather conditions, in the Q (7 mm) band using the standard VLA continuum mode 
(4IF, 50 MHz per IF). To minimize the effect of atmospheric phase noise, we observed in 
the fast-switching mode with a cycle of 120 seconds. The data were edited and calibrated 
by applying the complex gain solution from the calibration source following the standard 
VLA procedures. The flux density scale was determined by observing the source 3C286 for 
which we assumed flux densities of 1.45 Jy at 7 mm. In order to correct the amplitude 
and phase of the interferometer data for atmospheric and instrumental effects, we observed 
before and after every on-source scan the calibrator J1851+005. Standard calibration and 
data reduction were performed using AIPS. Maps were made by Fourier transformation of 
the interferometer data using the AIPS task IMAGR. The noise level in the images was in the 
range 0.12-0.16 mJy beam~^. The array was in the D configuration, resulting in a synthesized 
beam of ~ 2" at 43.4 GHz. In this observing session the source IRAS 18566+0408 was also 
observed, finding that its 7 mm emission probably arises from dust. This source is discussed 
in detail in Araya et al. (2007). 

We also made use of the archival VLA radio continuum data and found that all three 
regions have been previously observed at cm wavelengths. A summary of the archive data 
used here is given in Table [TJ 



^ The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 



2. 



Observations 



Astronomy Observatory on 2003 April 25. Each 




source was observed for 45 minutes, under 



-4- 



3. Results and discussion 

We detected 7 mm continuum emission towards each of the three observed sources. The 
position, flux densities, and angular size of the detected radio sources are given in Table [2j 
In what follows we present the results of our continuum observations and discuss the nature 
of the detected sources in each of the regions, individually. For the derivation of physical 
parameters we used the distances given by Sridharan et al. (2002) and listed in Table 3. 

3.1. IRAS 18470-0044 

Figure [1] shows a contour map of the 7 mm continuum emission observed towards IRAS 
18470-0044. The emission arises from two distinct compact sources separated by 28" (labeled 
A and B). Both sources were also detected at 3.6 cm and 6.0 cm wavelengths. A contour 
map of the 3.6 cm emission is shown in the lower panel of Fig. [H Component A is associated 
with a massive dust core detected at 850 and 450 fim by Williams et al. (2004) and at 1200 
yum by Beuther et al. (2002a). At 1.2 mm the core has major and minor axis of 17 and 
12", respectively, implying a core radius, estimated from the geometric mean of the axis, of 
0.29 pc (assuming a distance of 8.2 kpc). The mass of the core determined from the 850 /im 
observations is 720 Mq. The peak position of the core is marked with a cross in Fig. [TJ 

The spectral index of the emission between 5.0 and 43.4 GHz are — O.liO.l and — O.liO.l 
for the east (A) and west (B) components, respectively, indicating that the emission in this 
frequency range is optically thin thermal emission. We conclude that the emission from 
these two objects is free- free emission from ionized gas. Table[3]hsts the distance [col. (2)], 
the derived parameters of the regions of ionized gas: diameter [col. (3)], emission measure 
[col. (4)], and electron density [col. (5)], the minimum number of ionizing photons required 
to maintain the ionization of the nebula [col. (6)], and inferred spectral type of the exciting 
star [col. (7)]. They were calculated following the formulation of Mezger and Henderson 
(1967), assuming that the gas has constant electron density and an electron temperature of 
10^ K. 

The regions of ionized gas have radii of 0.021 and 0.029 pc, in the range of those of UC 
HII regions, and densities of 4.3 x 10^ and 2.5 x 10^ cm~^, much lower than those of UC 
HII regions. We suggest that these small regions of ionized gas are deeply embedded within 
molecular cloud cores with high densities and large turbulent motions, and have already 
reached pressure equilibrium with the dense ambient gas (e.g., De Free, Rodriguez, & Goss 
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1995; Xie et al. 1996). The equilibrium radius is given by (see Garay & Lizano 1999), 

where no is the density of the ambient gas, is the rate of ionizing photons emitted by the 
exciting star, and A?; is the hne width of the molecular emission from the ambient gas. The 
time to reach pressure equilibrium is 

.o .^J \'/V4x lO^cm-^V/V \'/V3kms-i\'/' 

Component A is projected right at the center of the massive and dense core which has an 
average molecular gas density of 1.4 x 10^ cm~^ and a line width in optically thin molecular 
transitions of 2.8 km s~^ (Beuther et al. 2002a). Using these values and the total number of 
ionizing photons per second needed to excite this component, of 2.1x10^^ s~^, we find that 
the equilibrium radius is ~ 0.014 pc, in good agreement with the observed radius of 0.021 
pc. The time to reach pressure equihbrium is ~ 2 x 10^ yrs. The actual age of the IRAS 
18470-0044 massive star forming region may be much larger than this value and thus might 
not correspond to an object in a very early stage of evolution. 

The total far- infrared luminosity of the IRAS 18470-0044 source, whose angular extent 
encompasses both regions of ionized gas, is 7.9 X 10^ Lq (Sridharan et al. 2002). On the 
other hand, the total luminosity of the ionizing sources is 2.3 x 10^ Lq. The difference in the 
inferred luminosities from the radio and FIR observations is most likely due to the presence of 
a group of non-ionizing stars which account for a fraction of the total bolometric luminosity 
but produces a negligible amount of ionizing photons. This explanation is strongly supported 
by our analysis of unpublished archive Spitzer IRAC observations which show the presence of 
several embedded sources within the angular extent of the IRAS source. However, it should 
be noted that the presence of dust absorbing ionizing stellar UV photons is most probably 
also contributing to the difference in the luminosities inferred from the radio and the FIR 
observations. 



3.2. IRAS 19217-1-1651 

IRAS 19217-1-1651 is associated with a massive dust core detected at 850 and 450 //m 
by WiUiams et al. (2004) and at 1200 /im by Beuther et al. (2002a). At 1.2mm the core has 
major and minor axis of 18 and 15", respectively, implying a core radius, estimated from the 
geometric mean of the axis, of 0.42 pc (assuming a distance of 10.5 kpc). The mass of the 
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core determined from the 850 /im observations is 2 x 10^ Mq. Beuther et al. (2002b) detected 
class II methanol and water masers near the peak position of the dust core, suggesting that 
this is a massive star forming region in an early stage of evolution. 

The upper panel of Fig. [2] shows a contour map of the 7 mm continuum emission 
observed with the VLA towards IRAS 19217+1651. Most of the emission arises from a 
single bright source, although there is evidence of a weaker component to the north. VLA 
observations at centimeter wavelengths (X, U, and K bands) show that emission arises from 
two components. This is illustrated in the bottom panel of Fig. [2] which presents a contour 
map of the emission observed at 22.5 GHz, showing that the emission arises from a bright 
compact component (labeled A), with an angular size of ~ 0.16", and a more extended, 
weaker component (labeled B) located ~ 1" NE of component A, with an angular size of 
0.8". The peak of the 7 mm source coincides with component A. The cross in both panels 
marks the peak position of the large dust core. 

The radio continuum spectra of component A, in the range from 8.4 to 43.4 GHz, is 
shown in the upper panel of Figure HI The spectra is reasonably well modeled (dotted line) 
by that of an homogeneous region of ionized gas with an emission measure of 7.0 x 10^ pc 
cm~^ and an angular size of 0.19". The region of ionized gas is optically thick below 15 
GHz. The total number of ionizing photons per second required to excite this region is 
5.5x10^^ s~^, which could be supplied by an 09.5 ZAMS star. The total luminosity of the 
region inferred from the radio observations is 5.4 x 10^ Lq, which is close to the luminosity 
inferred from the IRAS observations of 7.9 x 10^ Lq. Possible explanations for the difference 
in luminosities is that a fraction of the ionizing stellar UV photons is absorbed by dust 
within the ultracompact HII region or that non-ionizing stars are present in the region (or a 
combination of both). Our analysis of unpublished archive Spitzer observations show a single 
bright source associated with IRAS 19217+1651 that is detected in all four IRAC bands. 

The radio continuum spectra of component B, in the range of 8.4 to 22.5 GHz, is 
consistent with a total flux density of 6.6+1.0 mJy at 43.4 GHz, leaving a flux density 
of 43.8+1.0 mJy for component A at this frequency. The flat spectrum of component B 
is indicative of optically thin free-free emission. The direct interpretation of these data is 
that, since the total number of ionizing photons per second required to excite this region is 
6.9x10^^ s~^, we are observing an independent HII region ionized by a BO. 5 ZAMS star with 
a luminosity of 2.0 x 10"^ Lq. However, the high angular resolution images (see the 22.5 GHz 
image in Fig. 2) suggest that component B could be ionized gas outflowing from component 
A. Is this alternative possible? The angular separation between components A and B is ~ 1" 
that at a distance of 10.5 kpc gives a physical distance of 1.6 x 10^^ cm. Moving at a velocity 
of 10 km s~^, the ionized gas will take about 5,000 years in moving from component A to 
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component B. On the other hand, since component B has a deconvolved angular diameter 
of 0''8, we estimate for this source an average electron density of 8.9 x 10'^ cm~^ and a 
recombination timescale of only 14 years. Then, we conclude that component B requires its 
own ionizing star. 

3.3. IRAS 23151+5912 

IRAS 23151+5912 is associated with a massive dust core detected at 850 and 450 /im 
by Williams et al. (2004) and at 1200 /xm by Beuther et al. (2002a). At 1.2mm the core 
has major and minor axis of 16.5 and 14.4", respectively, implying a core radius of 0.21 
pc (assuming a distance of 5.7 kpc). The mass of the core determined from the 850 //m 
observations is 3 x 10^ Mq. Water masers were detected near the peak position of the core 
by Beuther et al. (2002b). 

The 7 mm continuum observations towards IRAS 23151+5912 shows the presence of a 
single source, as illustrated in the upper panel of Figure [31 We also detected emission from 
this source at 3.6 cm (with a total flux density of 0.27+0.04 mJy), as shown in the lower 
panel of Figure 3. The radio continuum spectrum, shown in the bottom panel of Fig. |H can 
be well fitted by that of an uniform density region of ionized gas with an emission measure 
of 2.3 X 10^ pc cm~^ and an angular size of 0.02"(or 6 x 10~^ pc). The high emission measure 
and small size of this HII region are characteristics of hypercompact HII regions (Kurtz 
2000, Hoare et al. 2007). Unfortunately, the modest angular resolution of our data does 
not allow a direct measurement of the angular size of this source and only an upper limit is 
possible (see Table 2). Additional observations of very high angular resolution are needed to 
test the nature of this source as a hypercompact HII region. Hypercompact HII regions are 
thought to mark the earliest stages of evolution of regions of ionized gas, being formed in the 
accretion phase of hot molecular cores (Keto 2002, 2003; Gonzalez-Aviles, Lizano & Raga 
2005; Avalos et al. 2006). We conclude that IRAS 23151+5912 is indeed associated with 
a massive star forming region in a very early stage of evolution. We note that the spectral 
index of the 7 mm source between 8.4 and 43.4 GHz is 1.1 +0.2. This value is typical of the 
spectral index of hypercompact HII regions at centimeter wavelengths, which suggests the 
presence of non-uniform density gas. With the presently available data we cannot however 
discern whether or not the ionized gas exhibits density gradients. In particular, we cannot 
rule out the possibility that the source is an ionized stellar wind. 

The rate of ionizing photons needed to excite the HII region is 5.5x10^^ s~^, which could 
be supplied by an BO. 5 ZAMS star with a luminosity of 1.1 x lO'* Lq. The FIR luminosity 
of the region inferred from the IRAS observations is, however, 1.0 x 10^ Lq (Sridharan et 
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al. 2002). The most likely explanation for the discrepancy in the luminosities inferred from 
the radio and FIR observations is that, due to the coarse angular resolution of the IRAS 
observations, the FIR luminosity includes the contribution of the star or stars ionizing the 
large cometary HII region located 30" NW of the dense core (see Fig. 5). This cometary HII 
region has a flux density of ~20 mJy, which requires a BO. 5 ZAMS star with a luminosity of 
2.0 X 10^ Lq to maintain its ionization. Another contributing effect could be absorption of 
UV photons by dust within the hypercompact HII region. There is no archive Spitzer data 
for this source. 



4. Summciry 

We made 7 mm, high-angular resolution, continuum observations, using the VLA, to- 
wards three regions of massive star formation thought to be in early stages of evolution. 
These regions are associated with energetic molecular outflows and two of them were not 
previously detected in radio continuum observations with the VLA. The main objective was 
to search for the presence of disks and/or hypercompact regions of ionized gas. Our main 
results and conclusions are as follows. 

Emission at 7 mm was detected toward the three regions. Towards IRAS 18470-0044 we 
detected two objects with spectral indices between 4.8 and 43.4 GHz of ~ —0.1, indicating 
that the emission is free- free radiation arising from optically thin HII regions excited by stars 
with spectral types BO. 5. We conclude that these small HII regions are embedded within 
massive and dense cores and have already reached pressure equilibrium with their dense and 
turbulent molecular surroundings. 

In the case of IRAS 19217-1-1651 we detect two components: component A is a hyper- 
compact HII region that is optically thick below 15 GHz, while component B is optically thin 
in the 8.4 to 22.5 GHz range. We discuss if component B could be ionized gas outflowing 
from component A but conclude that each component has its own ionizing star. 

Towards IRAS 23151-1-5912 we detected a hypercompact HII region that can be modeled 
as a source with a radius of 3 x 10~^ pc (or 57 AU) and an emission measure of 2.3 x 10^ 
pc cm~^. Wc conclude that IRAS 23151+5912 is indeed a massive star forming region in 
a very early stage of evolution, with the hypercompact HII region being formed while the 
suspected hot molecular core around it is still probably undergoing accretion. 

GG gratefully acknowledges support from the Chilean Ccntro dc Astrofisica FONDAP 
No. 15010003. LFR is thankful to the support of CONACyT, Mexico and DGAPA, UNAM. 
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Table 1. VLA ARCHIVE OBSERVATIONS ANALIZED 



IRAS source Project Observed date Preq. Conf. Synthesized 
name (GHz) Beam ("x") 



18470-0044 


AH810 21-MAR-2005 


4.86 B 


1.96 X 1.43 


1 O /( Tn f\f\A A 

lo4/U-UU44 


AS643 02- JUL- 1998 


8.46 BnA 


1.01 X 0.77 


19217+1651 


AS643 02- JUL- 1998 


8.46 BnA 


1.63 X 0.72 


19217+1651 


ABllOl 08-NOV-2003 


14.94 B 


0.40 X 0.36 


19217+1651 


ABllOl 08-NOV-2003 


22.46 B 


0.27 X 0.26 


23151+5912 


AR304 Ol-NOV-1993 


8.44 D 


8.98 X 6.82 


Table 2. 


OBSERVED PARAMETERS OP 7 mm EMISSION 


Radio source 


Peak position 


Plux density" 


Angular size 




a(2000) 5(2000) 


(mJy) 


(") 


18470-0044 A 


18 49 37.78 -00 41 01.8 


2.11+0.16 


1.2 X 0.9 


18470-0044 B 


18 49 35.94 -00 41 05.9 


2.45+0.14 


1.5 X 1.3 


19217+1651 


19 23 58.82 16 57 41.2 


50.4+0.16 


0.6 X 0.3 


23151+5912 


23 17 20.90 59 28 47.8 


1.70+0.19 


< 1 



^Plux densities are corrected for the primary beam response. The flux 
density of IRAS 19217+1651 is the sum from components A and B. 
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Table 3. DERIVED PARAMETERS 



Source 




D 


Diameter 


EM 


He 




S.T. 






(kpc) 


(pc) 


(pc cm""^) 


(cm~^) 






18470-0044- 


■A 


8.2 


0.042 


1.1 X W 


4.3 X 10^ 


2.1 X 10^6 


B0.5 


18470-0044- 


B 


8.2 


0.057 


5.3 X 10^ 


2.5 X 10=* 


1.5 X 10^*5 


B0.5 


19217+1651 


-A 


10.5 


0.009 


7.0 X 10^ 


2.5 X 10^ 


5.5 X 10^^ 


09.5 


19217+1651 


-B 


10.5 


0.041 


4.3 X 10^ 


8.4 X 10^ 


6.3 X 10^^ 


BO 


23151+5912 




5.7 


0.0006" 


2.3 X 10" 


1.7 X 10^ 


6.2 X 10^^ 


B0.5 



^The size of this source is estimated from the fit to its spectrum and not from 
an observed size. 
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Fig. 1. — VLA radio continuum maps towards IRAS 18470-0044. Beams are shown in the 
lower left corner of each panel. Top : 43.4 GHz map. Contour levels are -1, 1, 2, 3, 4 and 5 
X 0.25 mJy per beam. Beam is 2.1 x 1.6 ". Bottom : 8.46 GHz map. Contour levels are -1, 
1, 2, 3, 4, 5 and 6 x 0.20 mJy per beam. Beam is 1.0 x 0.77 The cross indicates the peak 
position of the dust core (WiUiams et al. 2004). 
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Fig. 2. — VLA radio continuum maps towards IRAS 19217+1651. Beams are shown in the 
lower left corner of each panel. Top : 43.4 GHz map. Contour levels are -1, 1, 2, 3, 5, 10, 
20, 30, 50 and 70 x 0.50 mJy per beam. Beam is 1.42 x 1.35 Bottom : 22.5 GHz map. 
Contour levels are -1, 1, 2, 3, 5, 10, 20, 40, 60, 90 and 120 x 0.30 mJy per beam. Beam is 
0.27 X 0.26 ". The cross indicates the peak position of the dust core (WiUiams et al. 2004). 
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Fig. 3. — VLA radio continuum maps towards IRAS 23151+5912. Beams are shown in the 
lower left corner of each panel. Top : 43.4 GHz map. Contour levels are -1, 1, 2, 3 and 4 x 

0. 40 mJy per beam. Beam is 1.94 x 1.52 Bottom : 8.44 GHz map. Gontour levels are -1, 

1, 2, 3, 4, 5, 7, 9, 11 and 13 x 0.060 mJy per beam. Beam is 9.0 x 6.8 The cross indicates 
the peak position of the dust core (Williams et al. 2004). 
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Fig. 4. — Radio continuum spectra of the regions of ionized gas with high emission measures. 
Dotted hues: fits to the observed spectra with a theoretical model of an homogeneous region 
of ionized gas. The fitted EM (in pc cm~^) is given in the upper right corner. Top: Spectra 
of radio component A towards IRAS 19217+1651. Bottom: Spectra of the radio source 
towards IRAS 23151+5912. 
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Fig. 5. — VLA radio continuum 8.44 GHz image towards IRAS 23151+5912. The beam is 
shown in the lower left corner. Contour levels are -1, 1, 2, 3, 4, 5, 7, 9, 11, 13, 15, 17, and 
21 X 0.060 mJy per beam. Beam is 9.0 x 6.8 ". The cross indicates the peak position of the 
dust core (WiUiams et al. 2004). 



